First Order Premelting Transition of Vortex Lattices 
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Vortex lattices in the high temperature superconductors undergo a first order phase transition 
which has thus far been regarded as melting from a solid to a liquid. We point out an alternative 
possibility of a two step process in which there is a first order transition from an ordinary vortex 
lattice to a soft vortex solid followed by another first order melting transition from the soft vortex 
solid to a vortex liquid. We focus on the first step. This premelting transition is induced by 
vacancy and interstitial vortex lines. We obtain good agreement with the experimental transition 
temperature versus field, latent heat, and magnetization jumps for YBCO and BSCCO. 
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Phase transitions involving vortex lattices in the high 
temperature superconductors is an area of active study 
Below a critical value of the magnetic field, 
vortex lattices in YBaaCuaOr-^ (YBCO) and 
BizSrzCaCuaOg (BSCCO) §-§ undergo a first order 
phase transition. This conclusion comes from latent heat 
measurements Q as well as jumps in the resistivity [^,||J9| 
and in the magnetization It has generally been 

assumed that this is a melting transition from a vortex 
solid to a vortex liquid. In this paper we suggest the 
possibility that the melting transition actually occurs in 
two steps as the temperature increases; the first step is a 
first order premelting transition from an ordinary vortex 
lattice to a soft solid with a small but finite shear modu- 
lus, and the second step is the first order melting of the 
soft solid into a vortex liquid. In this paper we focus on 
the first step. We present an analytic theory of a first 
order premelting transition in which the shear modulus 
jumps discontinuously. The transition is induced by in- 
terstitial and vacancy line defects in the vortex lattice 
which soften the shear modulus cge ■ We find good agree- 
ment with the experimental curve of transition temper- 
ature versus field, latent heat and magnetization jumps 
for YBCO and BSCCO. In the soft solid phase the super- 
conducting phase coherence along the field is destroyed 
by the wandering of the defect lines which become en- 
tangled in the vortices of the soft solid lattice JlO| , pl| . 
However, since wandering is energetically costly, the su- 
perconducting correlation length along the c-axis is long. 
Finally we speculate about the relation between our pro- 
posed two-step transition and the well known peak effect 

Let us describe our scenario for premelting. Our ap- 
proach follows that of Granato who showed that intersti- 
tial atoms soften the shear modulus of ordinary crystals 
and lead to a first order transition We start with 

a vortex lattice in a clean layered superconductor with 
a magnetic field H applied perpendicular to the layers 
along the c-axis. We consider the vortices to be corre- 
lated stacks of pancake vortices. We will assume that 



the transition is induced by topological defect lines, i.e., 
vacancies and interstitials. In a Delaunay triangulation 
JlSf a vacancy or an interstitial in a triangular lattice is 
topologically equivalent to a pair of bound dislocations 
Jl6t as well as to a twisted bond defect |l7). High tem- 
perature decoration experiments and Monte Carlo 
simulations have found such defects to be thermally 
excited. The introduction of these defects softens the 
elastic moduli. Since the energy to introduce intersti- 
tials and vacancies is proportional to the elastic moduli, 
softening makes it easier to introduce more defects. The 
softening also increases the vibrational entropy of the vor- 
tex lattice which leads to a premelting transition. The 
transition is driven by the increased vibrational entropy 
of the vortex lines of the lattice, and not by the entropy 
of the wandering of the defect lines. In fact Frey, Nel- 
son and Fisher |lCf| showed that a phase transition driven 
by the entropy of wandering flux lines occurs at a much 
higher magnetic field than what is observed experimen- 
tally. In the vicinity of the experimentally observed first 
order phase transition, wandering in the transverse direc- 
tion by more than a lattice spacing is energetically quite 
costly and therefore rare. (The energy scale is set by eoS 
Here s is the interplane spacing and Eq, the energy 
per unit length of a vortex, is given by Eo = {(j)o/4:TrXab)'^ 
where (f>o is the flux quantum and Xab is the penetration 
depth for currents in the ab plane.) 

Experimentally the resistivity jumps up at the tran- 
sition from zero to a flnite value as the temperature in- 
creases. This is consistent with our model since the soft 
solid will have a finite resistivity due to the motion of in- 
terstitial (and vacancy) lines. The barrier for the motion 
of interstitials is very small |l§ and is of order 10 ^Eq 
per unit length, where Eo — 2eo. The defect lines act 
like a liquid of lines existing in a soft solid host. Notice 
that if one tries to measure the shear modulus of such 
a system using resistivity measurements, only the defect 
lines would move relative to the pinned soft solid and one 
would deduce the shear modulus was zero [ |T8| , p^|jl^ ] . 

The first order transition is nucleated in a small re- 
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gion by a local rearrangement of existing line segments. 
Slightly above the premelting temperature Tp a vortex 
line can distort and make an interstitial and a vacancy 
line segment that locally create a soft solid. This is the 
analog of a liquid droplet which nucleates melting of a 
crystal. The role of the surface tension is played by the 
energy to connect the interstitial segment to the rest of 
the vortex line. This connection can be a Josephson vor- 
tex lying between planes or a series of small pancake vor- 
tex displacements spread over several layers. When the 
length £ of the interstitial and vacancy segments equals 
the critical length ic, the energy gained by premelting 
equals the energy cost of the connections. When £ > £c, 
it is energetically favorable for the defect segments grow 
to the length of the system. 

To study premelting we assume that we have a vor- 
tex lattice with interstitial and vacancy lines extending 
the length of the lattice. Our goal is to find the free en- 
ergy density as a function of the concentration n of defect 
lines. The free energy density is f = fo + fw + fvib + fwan 
where fo is the free energy density of a perfect lattice, 
is the work needed to introduce a straight interstitial or 
vacancy line into the lattice, fmb is the vibrational free 
energy density of the system, and /u,an is the free energy 
due to the wandering of the defect lines over distances 
large compared to the lattice spacing. We now examine 
these terms in detail. 

fo, the free energy density of a perfect rigid flux lattice, 
IS given by the London term 0,|o): 

fo = ^ In ( I , « B « H,2 



(1) 



where B is the spatially averaged magnetic induction, ^ab 
is the coherence length in the ab plane, and ry is 0.130519 
for a hexagonal lattice and 0.133311 for a square lattice 
I p^ . For B near Hc2, fo is given by the Abrikosov free 
energy [|ll 
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where the Ginzburg-Landau parameter k — Xab/Cab, and 
the Abrikosov parameter is 1.16 for a triangular lat- 
tice and 1.18 for a square lattice. 

To calculate fyn,, we follow Bulaevskii et al. j2^. We 
denote the displacement of the i^th vortex pancake in the 
nth plane from its equilibrium position by u(n, r^) where 
u = {ux, Uy) and the pancake position r — {tx, Vy). The 
Fourier transform u(k, q) — X^m/ ^{^^ ^i^) exp[i(k • -|- 
qn)]. k = (kx,ky) and q is the wavevector along the 
c-axis. fvib = —{ksT /V)\TiZmh where V is the volume 
and the vibrational partition function Zyu, is given by 
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where we have divided by the area tt^^j, of the normal 
core of a pancake . uji and it/ are the real and imagi- 
nary parts of u(k, q) and ie{x, y). The elastic free energy 
functional associated with these distortions is 



Tel = ^Wo^^Uj('7,k)ayM*(g,k) 



(4) 
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where i and je{x,2/}, the volume per pancake vortex is 
Vo = s(t)o/B, and s is the interplane spacing. The k sum 
is over a circular Brillouin zone = AnB/ipo- The ma- 
trix Uij is given by aij = cskikj + {ceek'^ + cnQ^)5ij 
where c_b, cee, and C44 are the bulk, shear, and tilt mod- 
uli, respectively, cb — cii — cge for a hexagonal lattice. 

— 2{\ — cos qs) / . Diagonalizing a^- leads to 2 eigen- 
values: Ai{kq) = ciifc^ -I- C44Q2 and At = CQ^k'^ + 044(3^, 
where A is the diagonal matrix, the subscript £ denotes 
longitudinal and t denotes transverse. Using this, we can 
integrate over m in (|^); the remaining sums over k and q 
are done numerically. At low fields {b = B/Hc2 < 0.25), 
the elastic moduli are given by 
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where Ac is the penetration depth for currents along the 
c-axis, 7 — Xc/Xab is the anisotropy, and C = 1. At 
high fields (6 > 0.5) [|l]j2|, cqq is altered by the fac- 
tor C « (1 - 0.5k-2)(1 - 6)2(1 - 0.585 + 0.2952) and 
the penetration depths in cn and C44 are replaced by 
A2 — A2/(1 — 5) where A denotes either Xat or Ac. 
In addition the last two terms of C44 are replaced by 
i3(/)o/(167r2A2). These replacements guarantee that the 
elastic moduli vanish at Hc2- For YBCO the temper- 
ature dependence of the penetration depths and coher- 
ence lengths are given by A(T) = A(0)(1 - {T/Tc))-^/^ 
II and U{T) = - (T/Tc))-i/2^ respectively. 

For BSCCO whose premelting field is two orders of mag- 
nitude below Hc2, A2(r) = A2(0)/(1 - [T/Tef ] and 
C\(7^)=Ca6(0)/(l-(T/Tc)4) @]. 

The free energy density /„, due to the energy cost of 
adding a vacancy or interstitial vortex line is difficult to 
calculate accurately pO[ ]. However, we can write down a 
plausible form for f^, by noting that a straight line de- 
fect parallel to the c-axis produces both shear and bulk 
(but not tilt) distortions of the vortex lattice. For ex- 
ample, if a defect at the origin produces a displacement 
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u that satisfies V • u = VoS{r)/s where S{r) is a two di- 
mensional delta function, then ^^(k) = ika/k"^ [ [To| . In- 
serting this in (||), we find that fw = (cge -I-cb)/2 where 
cb — X]k'^s('? ~ Ojk). Generalizing this to allow for a 
more complicated distortion and for a concentration n of 
line defects, we write |14 



fu 



dn{aiCeQ + q;2Cs) 
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where ai and a2 are dimensionless constants. We ex- 
pect the isotropic distortion to be small, i.e., 012 ^ 1, 
and the shear deformation to dominate, i.e., ai ^ a2- 
Integrating over n allows the elastic moduli to depend 
on defect concentration. We will assume that cb is in- 
dependent of n since we believe that the bulk modulus 
of the vortex solid is roughly the same as that of the 
soft soHd phase. To find cee{n) ||lj], we use its definition 
C66 = d^f/de^ where e is the shear strain. Assuming 
that Cb has negligible shear strain dependence, we find 
ceein) = ceeiO) + ai J^{d^CG(i{n)/de^)dn or 
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If we shear the lattice in the ab plane along rows sep- 
arated by a distance d, the shear modulus must be 
periodic in displacements equal to the lattice constant 
flo. We describe this with the simplest even periodic 
function: C(,e{u) = ce(,{u = 0) cos(27ru/ao) = C66(e 
0) cos(27r(i£/ao) where e — u/d. Then d'^cm{n) / de^ = 
—(3cQQ{n), where (3 = An'^d^ /a^^. Combining this with 
(0), we obtain CQQ{n) — C66(0) exp(— ai/3n). Thus the 
shear modulus softens exponentially with the defect con- 
centration n. This softening lowers the energy cost to 
introduce further defects, and increases the vibrational 
free energy f^b when C66('t-) is used in a^ . Substituting 
CQ(,{n) into our expression (p[) for yields 
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The last term we need to consider is fwam the free 
energy due to the wandering of the defect lines over dis- 
tances large compared to the lattice spacing. We can 
estimate /^jan with the following expression pO| 
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where = 3 for a triangular lattice (BSCCO) and 
= 4 for a square lattice (YBCO). £z can be thought 
of as the distance along the z-axis that it takes for the 
defect line to wander a transverse distance of one lat- 
tice spacing Cq. To go from one vacancy or interstitial 
site to the next, the defect line segment must jump over 
the barrier between the two positions. This gives £z a 
thermally activated form: £z ^ toG-xjpl—E /kBT), where 



io ~ ao(ei/es)"^^^ and E « ao(eieB)"^''^. ei is the line 
tension and is given by ei (eo/7^) lii(ao/^ab). Numeri- 
cal simulations [lO[ | indicate that the barrier height e_B is 
small and we use cb = 2.5 • 10~^eo- fwan itself is quite 
small compared to the other terms because of the high 
energy cost of vortex displacements. For example, in the 
soft solid phase at the transition fwan is about two orders 
of magnitude smaller than or fmt ■ Thus the transition 
is not driven by a proliferation of wandering defect lines 
because near the transition the high energy cost of vortex 
displacements is not sufficiently offset by the entropy of 



the meandering line |10|. 

Before we plot / versus n, we note that the difference 
between B and H is negligible for YBCO but can be a 
significant fraction of the premelting field Hp for BSCCO. 
To find the value of B to use in the Helmholtz free energy 
density /, we minimize the Gibbs free energy density G, 
i.e., dG/dB = where G = / - B • H/47r. Since the 
concentration dependence of B is negligible, we find B 
for n = for each value of H and T. Typical plots of 
A/ = f{n) — /(O) — fw + ^fvib versus n are shown in 
the inset of figure 1. The double well structure of A/ is 
characteristic of a first order phase transition. The equi- 
librium transition occurs when both minima have the 
same value of A/. We associate the minimum at n = 
with the vortex solid and the minimum at finite n with a 
soft vortex solid that has a small but finite shear modu- 
lus. The defect concentration at the transition is only a 
few percent. At higher concentrations A/ increases with 
increasing n because introducing defects costs compres- 
sional energy which is proportional to the bulk modulus. 
Thus defects do not proliferate. As an estimate of the 
softness at the transition, for n — 5%, C66(n) ~ 0.2c66(0) 
for BSCCO. The strain field £^p{k) produced by the de- 
fect determines whether the shear modulus is zero in 
the high temperature phase For dislocation loops, 

s'^Q{k) is singular as fc — s- 0, and the shear modulus is zero 



For vacancy and interstitial lines £Q^(fc) is 



-a/3 

at fc = 

finite, and hence the shear modulus is nonzero. 

In Figure |l] we fit the experimental first order transi- 
tion curves in the H — T plane using 2 adustable param- 
eters: ai and a2- As expected, ai ^ 02 and a2 <C 1 (see 
Figure |l|). The geometrical quantity f3 can have several 
values for a given lattice structure, depending on which 
planes are sheared. We choose (3 = tt^ tan^ (p where (f) is 
the angle between primitive vectors. Decoration experi- 
ments on BSCCO find a triangular lattice so we use 
(j) = 60°. For YBCO we choose — 44.1° which is very 
close to a square lattice which has 4> — 45°. Maki ||2^] has 
argued that the d-wave symmetry of the order parame- 
ter yields a square vortex lattice tilted by 45° from the 
a— axis. Experiments ||2^-[2^ on YBCO find (j) ranging 
from 36° to 45°. 

We can calculate the jump in magnetization AM at 
the transition using AM = —dAG/dH\T=Tp- The jump 
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in entropy As is given by As = —VodAG/dT\H^Hp 
where As is the entropy change per vortex per layer. 
The results are shown in Figure |[ We have checked 
that our results satisfy the Clausius-Clapeyron equation 
As = —{voAB/ATT)dHp/dT. We obtain good agreement 
with experiment well below Tc- Near Tc it is thought that 
the entropy jump is enhanced by microscopic degrees of 
freedom which are not included in our model. 

We can compare our results with the Lindemann cri- 
terion by calculating the mean square displacement < 
|up > at the transition using eq. (||): < |up >= 
-i2kBT/vo)J2ai.q9\nZ„b/dA{akq) where A is the di- 
agonal matrix similar to aij and a labels the 2 eigenval- 
ues. Defining the Lindemann ratio cl by cj^ =< |up > 
/al, we find that cl « 0.25 for YBCO at ffp = 5 T and 
that CL « 0.11 for BSCCO at Hp = 200 G. Here we have 
used the same values of the parameters that were used to 
fit the phase transition curves in Figure These values 
of Cl are consistent with previous values [p]J^j3l[| . 

Experiments have found little, if any, hysteresis 
This is consistent with our calculations. We can bound 
the hysteresis by noting the range of temperatures be- 
tween which the soft solid minimum appears and the solid 
minimum disappears. Typical values for the width of this 
temperature range are 300 mK for YBCO at H = 5T and 
1.3 K for BSCCO at H ^ 200 G. Another measure of the 
hysteresis can be found in the plots of A/ versus n. The 
barrier height Vb between the minima is low {VbVo ^ 30 
mK) which is consistent with minimal hysteresis. 

In going from the normal metallic phase to the vor- 
tex solid, two symmetries are broken: translational in- 
variance and gauge symmetry which produces the su- 
perconducting phase coherence along the magnetic field. 
In the soft solid phase, longitudinal superconductivity 
is destroyed by the wandering of the defect lines which 
become entangled with the soft solid vortices. (A vor- 
tex solid with entangled vortex lines has been termed a 
supersolid |]lo| , p| .) Even though line wandering is ener- 
getically costly and therefore rare, it does occur. As a re- 
sult, the correlation length along the c-axis will be quite 
long. This is consistent with measurements in YBCO of 
the c-axis resistivity which find that there is loss of vor- 
tex velocity correlations for samples thicker than 100 /im 
p2|-^ . For samples thicker than the longitudinal corre- 
lation length, the loss of longitudinal superconductivity 
coincides with the premelting transition [^5| . This agrees 
with experiments which indicate that the loss of super- 
conducting phase coherence along the c-axis coincides 
with the first order transition p2|-^. 

Since the soft solid is a lattice with a few percent of 
defect lines, the Fourier transform of the density-density 
correlation function should exhibit Bragg peaks. Rel- 
ative to the ordinary vortex solid, the intensity of these 
peaks would be slightly diminished by the defect lines, so 
it would be difficult to detect the transition via neutron 
scattering. In going from the soft solid to the normal 



metallic state, translational invariance is regained by a 
first order melting transition. Thus there are two tran- 
sitions: the premelting transition and the melting of the 
soft solid. Melting is observable in small angle neutron 
scattering experiments p6| which see a rapid decrease in 
the intensity of the Bragg spots. The region of the phase 
diagram where the soft solid exists may be quite narrow, 
of order a degree or less in temperature [|l^ . There is the 
intriguing possibility that our scenario of two transitions 
may be related to the peak effect in which the critical cur- 
rent as a function of temperature or field is observed to 
sharply increase below the melting transition [ p^ . This 
increase is believed to result from the enhanced pinning 
of flux lines due to the softening of the shear modulus cge 

To summarize we have discussed the possibility that a 
vortex lattice melts in two stages. First it undergoes a 
first order premelting transition into a soft solid followed 
by another first order phase transition into a liquid. The 
premelting transition is induced by vacancy and intersti- 
tial vortex lines that soften the shear modulus and en- 
hance the vibrational entropy. The entanglement of these 
defect lines with the vortex lines of the soft solid leads 
to the loss of longitudinal superconducting phase coher- 
ence. However, the correlation length corresponding to 
longitudinal superconductivity is quite long because line 
wandering is energetically costly and therefore rare. We 
obtain good agreement with the experimentally measured 
curve of transition temperature versus field, latent heat, 
and jumps in magnetization for BSCCO and YBCO. The 
Lindemann ratio cl is ~ 11% for BSCCO and ^ 25% for 
YBCO. The hysteresis is small. 
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FIG. 1. First order phase transition curves of magnetic 
field versus temperature, for YBCO and BSCCO. Parame- 
ters used for YBCO are ai = 2.55, a2 = 0.01485, (j) = U.r , 
\at{0) = 11861 @, s = 12A, UiO) = 15A, 7 = 5, and 
Tc = 92.74 K. Parameters used for BSCCO are Qi = 1.0, 
02 = 0.00705, (j) = 60°, Xab{0) = 2000l, s = 14l, 
Ca(,(0) = 30A, 7 = 200, and Tc = 90 K. For BSCCO we 
use the low field form of the elastic moduli from (jH]) and for 
YBCO we use the high field form. For fo we use (|l|) for 
BSCCO and (|) for YBCO. (For BSCCO we plot B vs. T 
because that is what ref. Q measured.) The experimental 
points for YBCO come from ref. ^ and those for BSCCO 
come from ref. [Q. Inset: Typical A/ versus n. 

FIG. 2. (a) and (b): Entropy jump As per vortex per layer 
versus Tp at the transition for YBCO and BSCCO. The exper- 
imental points for YBCO are from [^ and those for BSCCO 
are from Q. (c) and (d): Magnetization jump AM versus 
Tp at the first order phase transition for YBCO and BSCCO. 
The experimental points for YBCO are from [^] and those for 
BSCCO are from [Q]. For the theoretical points the values of 
the parameters are the same as in Figure 1 for all the curves. 



5 



76 



700 



600 - 



O 500 

o 
o 

^ 400 



CO 

CO 300 



O 
CD 



200 - 



100 - 







30 



Temperature [K] (YBCO) 

80 84 88 



92 



+ 



o 



o 



CNJ 



40 
30 
20 



O 



CO 
CO 

10 

o 

< -10 



+ +- 



"++- 



o 



+ 



1 1 1 1 1 














m 






T>T^ - 




1,1, 


m 

1 





■+++, o 



+ 



+ 



o Theory 
+ Experiment 



+ 



0.00 0.02 0.04 0.06 0.08 

n 



+ 



+ 



++- 



50 70 

Temperature [K] (BSCCO) 



10 



8 



6 CD 
C0_ 



4 

O 

o 



- 2 



90 







YBCO 



BSCCO 



1.0 



CD 
>^ 

B 

O 0.5 

I- 



CO 



0.0 
0.4 



(a) 



+ + + 



+ 



+ 
+ 



+ 



+ 



+ 



♦ + 

+ + + 

+ + 



(b) 



+ Experiment 
♦ Theory 



+ 
+ 
+ 




(c) 



(d) 



^ 0.3 h 
05 

CD 

" 0.2 h 



0.1 
0.0 



+ 



++ ♦ 

+ 



+ 

+ 
+ 
+ 



+ 



+ 



78 80 82 84 86 88 90 

T [Kelvin] 



30 



50 70 

T [Kelvin] 



90 



